Experimentally scalable protocol for identification of correctable codes 
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The task of finding a coiTectable encoding tlrat protects against some physical quantum process is in general 
hard. Two main obstacles are that an exponential number of experiments are needed to gain complete infor- 
mation about the quantum process, and known algorithmic methods for finding correctable encodings involve 
operations on exponentially large matrices. However, we show that in some cases it is possible to find such 
encodings with only partial information about the quantum process. Such useful partial information can be 
systematically extracted by averaging the channel under the action of a set of unitaries in a process known as 
twirling. In this paper we prove that correctable encodings for a twirled channel are also correctable for the orig- 
inal channel. We investigate the particular case of twirling over the set of Pauli operators and qubit permutations, 
and show that the resulting quantum operation can be characterised experimentally in a scalable manner. We 
also provide a postprocessing scheme for finding unitarily correctable codes for these twirled channels which 
does not involve exponentially large matrices. 



Introduction — The coherent experimental manipulation 
of quantum systems, and their application to various quan- 
tum information tasks, confronts significant limitations in the 
presence of noise, and in particular, decoherence. The discov- 
ery of quantum error correction codes enables methods for 
overcoming these limitations whenever the decoherence satis- 
fies various well-defined sets of conditions. Specifically, the 
appliciability of particular codes depends crucially on the de- 
tais of the physical noise model affecting a particular system. 
The standard approach for experimentally characterizing the 
full noise model and then assessing the usefulness of a given 
code are costly procedures, in particular, requiring a number 
of experiments yj, |2[] that grows exponentially with the num- 
ber of subsystems and the (classical) post-processing of ma- 
trices 10, whose dimension grows similarly exponentially. 
These limitations make the standard approach infeasible for 
the kinds of quantum information systems required in prac- 
tice. In this paper we describe a method that overcomes both 
of these issues. 

First, we demonstrate that partial information about the 
noise is sufficient to construct correctable encodings. This 
is motivated by the recent demonstration that valuable partial 
information about the noise can be extracted by symmetrising 
or "twirling" the quantum channel 0, S]. Twirling con- 
sists of averaging the action of a sequence of gates over some 
distribution of unitaries, in effect randomising aspects of the 
channel. This leads to scalable protocols which answer ques- 
tions of practical interest about the noise. In particular, it is 
possible to experimentally estimate the average gate fidelity 
between a quantum operation and the identity map |4(], as well 
as to obtain information about noise coiTelations in a quantum 
process H). 

We then focus on a particular type of twirl, the Pauli twirl, 
which has a number of special features that make finding 
noiseless and unitarily correctable encodings more tractable 
than the general case. First, the partial information obtained 



via a Pauli twirl can be accessed with only a polynomial num- 
ber of experiments, making such experimental characterisa- 
tion strictly scalable. Then, we discuss how this partial infor- 
mation lends itself to an algebraic algorithm for finding cor- 
rectable encodings. This new approach does not require the 
manipulation of exponentially large matrices. We also discuss 
how these encodings can be verified experimentally in an ef- 
ficient manner. 

Twirled Quantum Operations — Consider a trace preserv- 
ing linear quantum operation h{p) = J2i jlxlijPiPPj act- 
ing on n qubits, where the Pi j G CP„ are all Pauli opera- 
tors acting acting on n qubits. Trace-preservation requires 
that J2i jixlijPjPi ~ 1 (in particular, trx = 1), while 
Hermiticity-preservation requires [x]ij = [x]ji- Completely- 
positive (CP) operations have the aditional requirement that 
X > 0, which follows from the fact that CP maps can be writ- 
ten as a Choi-Kraus sum A{p) = 'J2k ^kpA], ifioll . This im- 
plies that for CP maps, Vi [xlu > 0, so the [xlu can be 
interpreted as probabiUties. 

Quantum operations require 0(2^") parameters to be fully 
described, and thus require an exponential number of experi- 
ments to be fully characterised lll|2|]. Due to this exponential 
cost, it is impractical to obtain a complete description about 
noise and decoherence acting on even a moderately large sys- 
tem of qubits. Useful partial information about the noise can 
be obtained by averaging the action of the quantum opera- 
tion under the composition U o AoU^ for unitary operations 
U{p) = UpU'^ randomly chosen according to some distri- 
bution 10, 13]. This averaging is known as a "twirl", and the 
averaged channel A(p) = j dp{U)U o K o (p) is known as 
the "twirled channel". The case where the distribution over 
unitaries is discrete is of particular interest. In that case, the 
twirled channel is given by A(p) = ^ Pi{Ui)Ui o KoU] (p), 
where Vr{Ui) is a probability distribution over the Ui. This 
leads us to our first result. 

Theorem : Any correctable code for a twirled channel A is 
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a correctable code for the original channel A up to an aditional 
unitary correction. 

Proof; Without loss of generality, consider A(p) = 
^j^. AfepAj. and a twirl with unitaries {Uj} where Ui ~ 1. 
Any unitary twirl is unitarily equivalent to a twirl that in- 
cludes the identity, and this unitary equivalence leads to the 
aditional unitary correction. A set of Choi-Kraus operators 
for A is then {Uj AkUj}. The existence of a correctable 
code for Ha under the action of A, with recovery operations 
^(p) — X^m RmpRln^ corrcsponds to a projector P into a 
subspace of the form TIa ®'Hb such that PR„iUj AkUjP = 
RmU]AkU,P&nARmU]AkUAp e is satisfied 

for all j,k,m la]. Since Ui ~ i, this implies PRmA^P = 
RmAkP md RmAk\p £ Ia^SIHb) for all fc, m. Thus 
is also correctable under the action of A. □ 

We now consider some specific twirls which are known to 
yield twirled channels which can be characterized efficiently 
via experiments. 

Pauli Twirling — It has been demonstrated 101 that twirling 
a channel A by the Pauli group T„ yields the effective channel 
A of the form 

In other words, the off-diagonal elements of x are eliminated. 
Channels of this form are known as Pauli channels. From here 
on A will denote the result of Pauli twirling A, unless stated 
otherwise. 

Pauli channels have a number of useful properties which fa- 
cilitate the search for some types of correctable codes. How- 
ever, the description of a general Pauli channel still requires an 
exponential number of parameters, and such parameters are 
not realistically accessible due to this exponential overhead. 
Instead, one can consider an additional twirl by the group n„ 
consisting of all qubit permutations. The channel An result- 
ing from a combination of a Tn twirl and a n„ twirl is what 
we call a permutation invariant Pauli (PIP) channel. Such 
channels can be written in terms of Choi-Kraus operators 

An(p) = 51^-11 "FIT II-^^-.-w4wP^-.>>w.iw, (2) 

where w = {wx,Wy, Wz) labels the number of the X, Y and 
Z Pauli operators, i/^ labels the Wx + Wy + Wz qubits over 
which Pw,iy„,i„ acts non-trivially, and iw labels which single 
qubit Pauli operator act on each of the qubits. The number 
Kn of different labels w needed to describe such a channel is 
Kn = ^n^ + + + 1, while for a fixed w there are 
/-iT^, = f , ) different i^^, and for fixed w and i/^, 

there are Kl^ = {""^^Zf'"') (""+"0 different i^. 

Two equivalent descriptions of the channel (with the same 
number of parameters) are the Choi-Kraus decomposition and 
the diagonal representation. The Choi-Kraus decomposition 
(|2]i can be rewritten as An(/o) = P^^^wip)^ where A/^ 



are the superoperators 

M^iP) = 1717 E 1^ E ^w,.„,i„pPw,.„4w (3) 

The are trace-preserving channels which apply each of 
the Pw,i/w.iw for a given w with the same probability. The 
AJP^ form a basis for PIP channels. 

Note that, given some n-fold tensor product of Pauli op- 
erators Pw,i/„,iw' we have An(-Pw,i/w,iw) = ^wPw,u„,i„ for 
some real constant Aw G [—1, 1], since Pauli operators ei- 
ther commute or anti-commute. Moreover An is self-dual, 
i.e. An = a fact that follows directly from the definition 
of Pauli channels. Since the {Pw.:y„,i^} form an operator ba- 
sis, it follows that the Aw are the eigenvalues of An, with high 
degeneracy, as they depend only on w. Thus, Pauli channels 
are diagonalisable and Hermitian. 

If we define \A)){{B\p = ,^A tr(i?V) then the diagonal 
representation of the channel An in terms of its eigenopera- 
tors is An(/o) = I]w=o ^^M^ip)^ where are the super- 
operators 

M^ip) = EE l^w,.„,iw»((^w..„,iwlP- (4) 

The diagonal description of An allows for straighforward 
description of the composition of channels. For example, sup- 
pose two Pauli channels A[|^ and A[^^ are be composed to 
yield A^^"* o Aj^'. This simply translates to the multiplication 
of the corresponding eigenvalues for each of the two channels 
because all Pauli channels commute. 

It is crucial to note that, because there is at most a polyno- 
mial number of different eigenvalues, they can be estimated 
efficiently by determining how a Pauli observable in a w class 
is scaled under the action of the twirled channel, as described 
in Ref. |5]. Thus, all the eigenvalues of a PIP channel can be 
estimated experimentally in a scalable manner 

Unitarily Correctable Subsystems — Let A and B be sub- 
systems of a Hilbert space H ~ Ha ®Hb® Hk- Given 
some channel A : B{H) — > BJH), we say that Ha is a uni- 
tarily correctable subsystem ||9[] if there is a unitary recov- 
ery operation TZ acting on B{H) such that Vcta e B{Ha) 
and V(T_B S B{Hb) there exists a tb S B{Hb) such that 
TZ o K{ijA ® cfb) = cta ® tb- That is. Ha defines a cor- 
rectable encoding for A that can be returned to its initial loca- 
tion within the system Hilbert space with a single unitary op- 
eration. The problem of finding unitarily correctable subsys- 
tems for a unital channel is equivalent to finding the structure 
of the commutant of the noise algebra of A^ o A |@]. In terms 
of the Choi-Kraus operators {A^} for A, this "noise commu- 
tant" is defined as the set of operators that commute with the 
operators If A is unital and trace preserving, this 

commutant coincides with the fixed point set of A^ o A (the 
set of operators which are invariant under the action of this 
composed channel) jst]. This result can be refined somewhat 
for diagonalisable channels. 
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Proposition 1 — Let A be a unital, diagonalisable and trace 
preserving channel with eigenvalues and eigenoperators 
Li. Then, the noise commutant of At o A is the space spanned 
by eigenoperators Li with eigenvalues |Ai| = 1. 

Proof: A is unital and trace preserving, thus so is At o A. 
Moreover, A is diagonalisable, so At o A has eigenoperators 
Li with eigenvalues A* A; = |Aip, where the Ai are the eigen- 
values of A. Since At o A is unital, its fixed point set and its 
noise commutant coincide, and both are given by the space 
spanned by the eigenoperators L; with |Aip = 1. □ 

Pauli channels are unital channels, and since they are diag- 
onalisable and Hermitian, these channels have a particularly 
simple fixed-point set structure. In particular, we immedi- 
ately obtain the following result. 

Corollary 1 — Let A be a Pauli channel. Then the noise 
commutant of At o A is the space spanned by the eigenopera- 
tors with eigenvalues ±1. 

Before continuing, we discuss a special class of UCS codes 
considered in ll ill . We say that a UCS code is unitarily noise- 
less (UNS) for A if it is a UCS of A" for all n > 1, where 
A" is the channel A composed with itself n times. They are 
codes for which the recovery operation has the special form 
U = Ua ^TZb, that is, a unitary acting only on the subsystem 
in which information is preserved, and an arbitrary quantum 
channel on subsystem B. Interestingly, the sets of UCS and 
UNS codes coincide for Pauli channels. 

Corollary 2 — A UCS code of a Pauli channel A is also a 
UNS. 

Proof: As A is Hermitian, it commutes with its dual At = 
A, and thus we have that (A")^ o A" = (At o A)". Since the 
eigenvalues of A are real, this implies that the fixed point set 
of At o A is identical to the fixed point set of (A") ^ o A". □ 

PIP channel parameter space — Considering the Liouville 
representation of the superoperators {A/p,} and {M^}, it is 
easy to show that there is a linear invertible map : R^" 
R^" mapping the to the Aw ■ More explicitly, we have that 



Aw = Y.in\ 



(5) 



where (•, •) is the Hilbert-Schmidt inner product in the Li- 
ouville representation. This follows directly from the fact 
that (M^,M^) = (5w,v(Mw>^^w)- Similarly, because 
{MP, MP) = d^.^{MP,MP), we have 



{MP,M^) 
{Ml, Ml)' 



(6) 



The {pw} form a Kn — 1 simplex (the probability simplex, 
which is in fact the standard Kn — 1 simplex), and the {Aw} 
also form a Kn — 1 simplex, which we call the eigenvalue 
simplex. 

Explicitly computing the entries of the matrix from the 
matrix elements of Liouville representation of the superoper- 
ators is inefficient, as these representations are exponentially 
large in n. However, one can show that [il]v,w = 1— 2 ^^^,^^-' , 



where A'^(v, w) is the number of Pauli operators in the Kraus 
decomposition of an extremal channel (i.e. pv = 1 for 
some v) which anticommute with some Pauli operator corre- 
sponding to the equivalence class w [5]. General expressions 
for (v, w) can be obtained by simple counting arguments, 
which can be computed efficiently. 

The fixed points of An can be thought of as the observables 
that are conserved under the action of An in the Heisenberg 
picture, as An is self-dual. Once the fixed points of these 
channels are determined, in order to determine possible en- 
codings one needs to compute the possible algebras generated 
by these fixed points. Note that, because of the degeneracy 
of the eigenvalues of a PIP channel, the existence of a single 
weight class with eigenvalue 1 corresponds to a large number 
of Pauli operators which are fixed points of the channel. We 
now describe how this allows us to find correctable codes for 
a An. 

Finding Correctable Codes — The identification of the 
fixed point set of a unital channel can be used to find a UCS 
using the general algorithm described in Ref. |@, @]. How- 
ever, this algorithm requires the manipulation of exponentially 
large matrices, a problem shared with numerical algorithms 
used to search for more general codes Given the Pauli 
operators form an eigenbasis for PIP channels, the task of 
computing the algebra of conserved observables is relatively 
simpler than in the general case. One simply has to group the 
conserved Pauli observables into triplets of observables which 
satisfy the commutation relations for su{2). These commuta- 
tion relations can be computed without writing the observ- 
ables explicitly in a particular representation. By choosing 
the largest set S of mutually exclusive triplets which com- 
mute with each other, one implicitly describes how to encode 
a noiseless Hilbert space of dimension 2l'^l . Before discussing 
this algorithm, we point out the following result, which can be 
obtained by direct computation. 

Proposition 2 — Given PauU operators {Pj, Pk, Pi} sat- 
isfying the commutation relations [Pj,Pk] = ^jklPl^ 



where 



£jki-^, ^jki is the Levi-Civita symbol and 



^jki - ^]ki—, £jki 
Sj,Sk,si e {±l,±i}, then {sjPj,SkPk,siPi} obey the 
su{2) commutation relations. 

It is clear that {sjPj,SkPk,siPi} are unitarily related to 
{X, Y, Z}, as all Pauli operators have the same spectrum and 
Proposition 2 guarantees they have the same commutation re- 
lations. The unitary which performs the encoding is guaran- 
teed to be in the Clifford group, since it maps a set of Pauli 
operators to another set of Pauli operators with the same com- 
mutation relations. Standard techniques can be applied to de- 
termine which Clifford group operations implement a desired 
encoding 11211 . 

This leads to the following algorithm for finding a cor- 
rectable encoding given the eigenvalues of the PIP channel A: 
(i) Enumerate the Pauli operators with eigenvalues 1 under the 
action of A, and call thise set F. (ii) Choose a triplet of Pauli 
operators satisfying the commutation relations in Proposition 
2 - if none can be found, the search is over, (iii) Remove this 
triplet from F, as well as all operators that do not commute 
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with the triplet, and go back to step ii. The number of mutu- 
ally exclusive triplets found in this manner corresponds to an 
allowable number of encoded qubits which can be protected 
from the action of An- 

Finding unitarily correctable subsystems is similarly sim- 
ple. One can easily compute the conserved observables of the 
channel Ajj o An 10 ■ In the case of PIP channels, this cor- 
responds to finding observables with eigenvalue ±1 for the 
channels An, so that these observables have eigenvalue 1 for 
A\j o An = A^. The same procedure described above can be 
applied to the ±1 eigenspace to find a correctable encoding. 

Example 1 — Consider the 2 qubit PIP channel with Kraus 
operators proportial to {11, ZZ^. The Pauli operators with 
eigenvalue 1 are It, XX,YY,ZZ,XY,YX, and 1Z,Z1. 
Out of this set, {XX, XY, IZ} satisfy the commutation rela- 
tions, and no other triplets which commute with these can be 
found, so a single qubit can be encoded noiselessly through 
this channel. □ 

Example 2 — Consider the 2 qubit PIP channel with Kraus 
operators proportional to {li,YX,XY}. The eigenopera- 
tors with eigenvalue 1 are 11, XY, YX, and ZZ. There are 
no triplets with the right commutation relations. The eigen- 
operators with eigenvalues —1 are IZ, Zl,XX,YY. If we 
consider the ± 1 eigenspace, we obtain the same eigenopera- 
tors with eigenvalue 1 as the previous example, and thus there 
exists a UCS consisting of a single qubit. □ 

In the case of the previous examples, we want to map the 
generating set of Pauli operators {iX, iZ} to the generating 
set {XX, tZ}, which can be done by a controlled-NOT gate, 
where the second qubit is the control, and the first qubit is the 
target. 

The question as to whether this algorithm can find all uni- 
tarily correctable encodings for a PIP channels remains open. 

Verification of UCS — The conditions for error correction 
are known to be robust against perturbations ifisll . and bounds 
on the fideUty of the encoded states can be computed given 
complete information about A. As argued earlier, this cannot 
be done in a scalable manner Moreover, our theorems apply 
to CP maps, while it is known that there are physical maps 
which are non-CP 11411 . For these reasons, it is important to 
check what the fidelity of the encoding constructed for An 
is under the action of A. This can be done by twirling the en- 
coded channel, and estimating its average gate fidelity with the 
identity channel. The estimate can be done efficiently, using a 
number of experiments which is logarithmic in the number of 
encoded qubits, as described in Ref. jst]. 

Discussion — In this paper we have shown that correctable 
encodings for a quantum operation can be found by search- 
ing for correctable encodings using the twirled version of that 



quantum operation. We investigated in detail the case of chan- 
nels twirled by Pauli operators and qubit permutations, and 
demonstrated a simple scheme for identifying encodings with 
unitary recovery operations. Such twirled channels are im- 
portant because they are described by a polynomial number 
of parameters which are experimentally accessible via a scal- 
able protocol. The scheme does not require the manipulation 
of exponentially large matrices, and the performance of the 
constructed correctable code can be estimated experimentally 
in an efficient manner Further work is needed to determine 
whether all unitarily correctable codes for PIP channels can be 
found through the scheme we described. Moreover, it would 
be interesting to investigate other twirls that yield more in- 
formation about the channel which may be used to find more 
correctable codes. 
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